Electronic states in 1/1 Cd 6 Yb and 1/1 Cd 6 Ca: Relativistic, correlation, and structural effects 
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The electronic structure of the rational approximants 1/1 CdeYb and 1/1 CdeCa to the stable icosahedral 
CdYb and CdCa quasicrystals is studied by the full-potential linear augmented plane wave method. A com- 
parison is made between several structural models. We show that the (relativistic) spin-orbit (SO) interaction 
and electronic correlations that are not described by the usual local density approximation, are essential for an 
accurate description of the electronic structure. In particular, we show that the SO interaction is responsible 
for a splitting of the Cd-4d and Yb-4/ peaks, and that the experimental peak positions can be reproduced by 
including a Hubbard U term in the Hamiltonian [U(Cd) = 5.6 eV, U(Yb) = 3.1 eV]. Our results show very 
good agreement with a photo-emission (PE) spectrum of 1/1 CdeYb [R. Tamura, Y. Murao, S. Takeuchi, T. Kiss, 
T. Yokoya, and S. Shin, Phys. Rev. B 65, 224207 (2002)] and a 350 eV PE spectrum of 1/1 Cd 6 Ca, which we 
present in this paper. Without the relativistic and correlation effects even a qualitative agreement with the PE 
spectra cannot be achieved. 

PACS numbers: 71.23.Ft,79.60.-i 



The recently discovered binary icosahedral quasicrystals 
(QCs) /-CdYb 1 and /-CdCa 2 have received considerable at- 
tention. Their electronic structures have been calculated and 
analyzed 3 from the viewpoint of their stability and differences 
with respect to ternary Al-based QCs. The purpose of this 
communication is to point out that two elements ignored in 
previously published calculations, namely, the spin-orbit (SO) 
coupling and correlation effects, are of crucial importance for 
an accurate description of the electronic structure of these 
QCs. In addition, this communication addresses the issue of 
the stability of various structural models proposed for the 1/1 
approximants of these QCs. 

The binary icosahedral QCs /-CdYb and /-CdCa were 
discovered 1,2 in 2000. Prior to this discovery the most studied 
QCs were the ternary alloys /-AlCuFe, /-AlPdMn, /-AlPdRe, 
d-AlNiCo, and J-AlCuCo. It has been anticipated 1 that binary 
QCs simplify the analysis of experimental properties. The 
study of /-CdYb and /-CdCa has indeed shed light on prop- 
erties that are arguably still not well understood in the above- 
mentioned ternary QCs, such as the negative temperature co- 
efficient of the resistivity 4 and the origin of the pseudogap in 
the electronic density of states (DOS)^ 

Theoretically the binary QCs require a more sophisticated 
approach than the above-mentioned ternary QCs because of 
the relatively heavy elements they contain: As we will show, 
the (relativistic) SO interaction and electron-electron interac- 
tions that are not described by the usual local density approx- 
imation (LDA), cannot be ignored. Because these effects are 
most prominent in /-CdYb, we will focus on this QC. 

1/1 CdeYb 6 is a rational approximant 7 to /-CdYb. Its struc- 
ture can be described 8 as a bcc packing of interpenetrating 
icosahedral clusters (Fig. [I]). /-CdYb is composed of the 
same clusters. 9,10 Whereas the exact structure of /-CdYb has 
not been determined, the structure of 1/1 CdeYb is known. 
Some uncertainty exists only about the positions of the four 
Cd atoms inside the so-called dodecahedral cavity (Fig. [TJ. 



1/1 CdeCa is isostructural with 1/1 CdeYb: Ca atoms replace 
Yb. In this paper we present an ab initio electronic structure 
study of 1/1 CdeYb and 1/1 CdeCa. Because of the structural 
proximity of the binary icosahedral QCs and their 1/1 approx- 
imants we expect that our results that do not depend on the 
long-range order of these materials, are also valid for /-CdYb 
and /-CdCa. 

We studied 1/1 CdeYb and 1/1 CdeCa with the computer 
program wlEN2k, n an implementation of the full-potential 
linear augmented plane- wave method 12 based on density func- 
tional theory. 13 For the exchange and correlation potential we 
used the LDA. 14 In our DOS calculations we included the SO 
coupling using the so-called "second variational treatment." 15 

Although the LDA gives reliable optimal atomic positions 
for many materials, for the same materials it may fail to pre- 
dict the proper electronic ground state. 16 In particular, the 
LDA description of the electronic structure of compounds 
containing (strongly correlated) transition-metal d or rare- 
earth/ electrons is often inadequate. 1 — As we will show be- 
low, a good example of this shortcoming of the LDA is the 
wrong prediction 3 for the position of the Yb-4/ states in the 
DOS of 1/1 CdeYb. A qualitative improvement is achieved by 
the LDA + U method, 18 which takes into account the orbital 
dependence of the Coulomb and exchange interactions, which 
is absent in the LDA J 7 - 

We applied the LDA + U formalism of Ref.[Hto the Cd-4d 
and the Yb-4/ states. The orbital potentials felt by these states 
depend on two adjustable parameters: the Coulomb parame- 
ter U (also called Hubbard U) and the exchange parameter J. 
For fully occupied orbitals like the Cd-4d and the Yb-4/ states 
in 1/1 CdeYb it is straightforward to show that the effect of U 
and J is to shift the orbital potential by —\{U — J). There- 
fore we could choose J = without loss of generality. We 
adjusted U to find the best possible agreement with a photo- 
emission (PE) spectrum of Ref.ll9land a PE spectrum that we 
present in this paper. Our optimal values are U (Cd) =5.6 and 



FIG. 1: (Color online) Unit cell of 1/1 Cd 6 Yb (spacegroup 123). The 
Cd (small) and Yb atoms (large) form a shell structure. In the cen- 
ter Cd atoms are located on four vertices of a cube. The next shell 
(the dodecahedral cavity) consists of Cd atoms placed on the ver- 
tices of a dodecahedron. The third shell is a Cd icosidodecahedron 
composed of triangles and pentagons. Yb atoms are located at the 
centers of the pentagons. The remaining Cd atoms and Cd atoms of 
the icosidodecahedra in neighboring unit cells (not shown) form a 
defect triacontahedron (Ref. 8). 



U(Yb) = 3.1 eV 

Different structural models have been proposed for 1/1 
CdeYb. They differ essentially only in the positions of the 
four Cd atoms inside the dodecahedral cavity (Fig. [I}. Struc- 
tural details are often important for the electronic structure. 
Therefore we made a comparison of the various structural 
models. In the original model 6 of 1/1 CdeYb the Yb and Cd 
atoms are on the same positions (scaled by the lattice parame- 
ter a) as Ru and Be in Ru3Bei7 (Ref. 20). Four additional Cd 
atoms are located 6 on four vertices of a cube [Fig. Eta)] in- 
side the dodecahedral cavity. The spacegroup of this original 
model is 123 (Ref. 0). 

To calculate a good electronic DOS it is in general neces- 
sary to use accurate atomic positions. 21,22,23 In particular, ide- 
alized positions may lead to a DOS with unphysical features 
(spikes), 21 a pseudogap at the Fermi energy Ep that is too 
narrow, 21 or shifted peaks. 22 Therefore we have performed 23 
a structural relaxation of the atomic positions of the original 
model of 1/1 Cd6 Yb. Our final positions- 3 - are relatively close 
to the original positions. 6 During the structural relaxation the 
total energy decreased by 17 meV / atom. In the DOS the 
peak around Ep — 9 eV became 14 % narrower [Figs. |3ja) 
and [3b)], indicating a reduced (spatial) overlap of the Cd-4d 
states. 

Use of the generalized gradient approximation 24 (instead 
of the LDA) resulted in marginally different atomic positions. 
The main changes in the DOS were a lowering of the Cd-4d 



FIG. 2: (Color online) Models of 1/1 Cd 6 Yb differ from each other 
in the positions of the Cd atoms inside the dodecahedral cavity (Fig. 
QJ. These positions and the unit cells are schematically indicated for 
models with symmetries (a) 123, (b) Pn3, (c) 1222, (d) C2, and (e) 
Imm2 (See text). 



states by 0.06 eV and an increase of the Yb-4/ peak by 0.01 
eV. In the present study we are not interested in such small 
differences. Therefore, we present results for the LDA only. 

Electrical resistivity, specific heat, and electron diffrac- 
tion experiments have suggested 25 that at low temperatures 
(T < 110 K) a structure with a simple cubic unit cell (space- 
group Pn3) is favored over the bcc 123 structure. (Later 
experiments 26 have indicated that the low-temperature phase 
is noncubic, but for computational reasons we will consider 
noncubic structures only in the body-centered models.) In 
the Pn3 model the Cd atoms inside the dodecahedral cavi- 
ties occupy vertices of cubes with an alternating orientation 
[Fig.Etb)]. The DOSs of the original 123 and Pn3 models are 
not significantly different. 27 We relaxed the Pn3 model. The 
total energy of the relaxed Pn3 model was 0.4 meV / atom 
lower than that of the relaxed 123 model. Further convergence 
studies would be necessary to establish whether this number 
is significant. For the present study, however, it suffices to say 
that we found that the DOSs of the relaxed 123 and Pn3 mod- 
els were practically identical, and that we therefore restricted 
ourselves to the computationally less involved body centered 
models, such as the 123 model [Fig.|3a)]. 

Another modification of the 123 model has been proposed 
on the basis of a recent x-ray diffraction study. 8 According to 
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FIG. 3: Electronic DOS of 1/1 Cd 6 Yb (spacegroup 123). (a) Original 
model (Ref. 6). (b)-(d) Relaxed model (Ref. 23). In (c) and (d) 
the SO interaction is included, (d) was calculated with the LDA + 
U method. All DOSs were convoluted with a Gaussian with a full 
width at half maximum of 30 meV. The DOS in (d) is not as smooth 
as the other DOSs, because fewer k points were used to calculate it. 



this study the central Cd atoms are located on vertices of a 
truncated cube [Fig.|2jd)]. When we tried to occupy four of 
these vertices, we found that there are only two inequivalent 
ways to do this (we discarded Cd-Cd nearest neighbor dis- 
tances that were shorter than 2.4 A). The first possibility is that 
the Cd atoms are on vertices of a rectangular parallelepiped 
[Fig.|2jc)]. This structure has spacegroup 1222. When we re- 
laxed these atomic positions we found that the atoms moved 
to the positions corresponding to the 123 symmetry [Fig.|2ja)], 
indicating that the 1222 structure is not stable. The second 
possibility is indicated in Fig. EJd). The spacegroup is C2. 
When we relaxed the atomic coordinates of this model we 
found that the Cd atoms moved to positions somewhere in be- 
tween vertices of a cube octahedron [Fig.|3e)] and vertices of 
an icosahedron. The spacegroup of this model is Imm2. The 
Imm2 model had a total energy that was 2 meV / atom lower 
than that of the relaxed 123 model, indicating that the Imm2 
structure is more stable than the 123 structure. The peak posi- 
tions in the DOSs of the relaxed Imm2 and 123 models were 
equal within 0.01 eV. 

Summarizing our results so far, we found that the different 
models of 1/1 CdeYb had similar electronic DOSs. The only 
significant change was a narrowing of the Cd-4d peak after the 
structural relaxation of the original 123 model. It is of particu- 
lar interest that the main features in the DOS were found to be 
insensitive to the precise positions of the four Cd atoms inside 
the dodecahedral cavity, which is consistent with the theoret- 
ical result of Ref. that these Cd atoms are weakly bound 
and do not play an important role in stabilizing 1/1 Cd6 Yb. In 
the rest of this paper we present results for the model with the 
highest symmetry, i.e., the 123 model. 

We calculated the electronic DOSs of 1/1 Cd 6 Yb and 1/1 
CdeCa. A separate structural relaxation was performed for 1/1 
CdeCa. Figured shows our results for 1/1 CdeYb obtained at 
various levels of sophistication. The result of Fig. Eta) was 
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FIG. 4: (Color online) Comparison of the calculated DOS of 1/1 
CdeYb with a PE spectrum (Ref. 19). The DOS was convoluted with 
a Gaussian with a full width at half maximum of 200 meV. 



obtained without SO coupling or orbital potentials (Hubbard 
U s) for the original 123 model. It is essentially equal to the 
result of Ref.|3 The main features are the Cd-4d peak around 
Ep — 9 eV, one Yb-4/ peak just below the Fermi energy, and a 
pseudogap near Ep. In Ref. 3 it has been shown that this sta- 
bilizing pseudogap is formed as a consequence of hybridiza- 
tion between the free-electronlike s and p states and low-lying 
(unoccupied) Yb-5d states. 

The DOSs of Ref. 3 and Fig.Et a ) are in disagreement with 
a PE spectrum 19 of 1/1 CdeYb, in which there are two Yb-4/ 
peaks around 2.0 and 0.7 eV below the Fermi energy. The 1.3 
eV split between these peaks is attributed 19 to the SO interac- 
tion. It has been suggested 19 that differences between the Yb- 
4f peak positions in the DOS of Ref .0 and in the PE spectrum 
of Ref. [jjj are due to a shape approximation for the potential 
that has been made in Ref. U The fact that our full-potential 
calculation is essentially equal to the result of Ref. yj shows 
that this is wrong. In fact, a comparison of Figs. Ob) and0c) 
shows that the most obvious discrepancy between Fig. Et a ) 
and the experiment, 19 namely, the absence of the SO splitting 
in Fig.[3ja), is a result of the absence of the SO interaction in 
the calculation. 

As mentioned above the structural relaxation of the 123 
model of 1/1 CdeYb led to a narrowing of the Cd-4d peak 
[Figs. Eta) andEtb)]. The SO interaction split the Yb-4/ and 
Cd-4d peaks [Fig. Etc)]. By including a Hubbard-/7 term in 
the Hamiltonian for the Cd-4d and Yb-4/ states we could shift 
these peaks to lower energies [Fig.EJd)]. We adjusted the U s 
to get the best possible fit with the PE spectra shown in Figs. 
|4|and|5| Our final DOSs reproduced the experimental peak 
positions very well (Figs.|4|and|51). 

The PE spectrum of 1/1 CdeCa (Fig. EJ) was measured 
by x-ray PE spectroscopy with an incident photon energy of 
350 eV at the temperature of 20 K using a GAMMADATA- 
SICENTA SES2002 hemispherical analyzer. Details of the 
sample preparation have been given elsewhere. 5 The Cd-4d 
states are greatly enhanced in the experimental spectrum since 
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FIG. 5: (Color online) Comparison of the calculated DOS of 1/1 
CdeCa with a PE spectrum (See text). The DOS was convoluted 
with a Gaussian with a full width at half maximum of 400 meV. 



TABLE I: Calculated and measured electronic contributions to the 
specific heat. The values are given per "mole formula unit", viz., 
CdeYb and Cds.rYb for the 1/1 approximant and the QC, respec- 
tively. 

7(mJmor 1 K -2 ) 
1/1 Cd 6 Yb /-CdYb 



Calculation / Experiment 



LDA calculation [Fig. Ha)] 

LDA calculation [Fig. 0b)] 

LDA + SO calculation [Fig. 0c)] 

LDA + SO + U calculation [Fig.0d)] 

Experiment (Ref.^fr) 

Experiment (Ref. 29) 

Experiment (Ref. 30) 



7.3 
7.4 
11.8 
6.1 



51 



19.2 
7 

7.5 



at this photon energy the cross section of the Cd-4d states 
is the largest among all the states that may contribute in this 
binding energy range: For instance, the cross section of Cd-4d 
is one or two orders of magnitude larger than those of Cd-5s 
and Ca-4s. 

The SO interaction and the strong electronic correlations 
in 1/1 CdeYb are not only important for the peak structure 



in the DOS, they are also essential for a proper description 
of the DOS near the Fermi energy. This is illustrated by a 
comparison of the electronic contributions to the specific heat, 
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(1) 



calculated at various levels of sophistication (Tabled, where 
we assumed that the phonon-correction factor A = 0. (For a 
proper comparison with experiment A can, of course, not be 
ignored. 28 ) Experimental values of 7 are also given (Table|I|>. 
It is remarkable 30 that 7 is much lower in the QC than in the 
1/1 approximant (Table This indicates a substantial rear- 
rangement of states leading to a considerable lowering^ of 
the DOS at the Fermi energy in /-CdYb. This rearrangement 
of states could be due purely to geometrical effects (Hume- 
Rothery), or geometrical effects coupled with increased hy- 
bridization between the Yb-5d states and the free-electronlike 
s and p states, a scenario that still remains to be explored. 

In 1/1 Cd 6 Yb the Yb-4/ states are fully 28 occupied. 
The states at the Fermi energy have essentially no Yb-4/ 
character. 28 Therefore, the Yb-4/ states do not contribute 2 ^ to 
transport properties or to the temperature dependence of ther- 
modynamic properties. This is in clear contrast to the mixed- 
valent compound YbCuAl, in which the Yb-4/ electrons are 
important for the heavy fermion character^! 

In conclusion, the four Cd atoms inside the dodecahedral 
cavity occupied vertices of a (deformed) cube octahedron in 
the most stable model of 1/1 CdeYb. The main features in the 
electronic DOS were insensitive to the model used. To get a 
good agreement with PE spectra it was necessary to include 
the SO interaction in the Hamiltonian, which is responsible 
for the splitting of the Cd-4d and the Yb-4/ peaks in the DOS. 
In addition, the Cd-4d and Yb-4/ electrons have strong corre- 
lations that cannot be treated in the LDA. We could describe 
these correlations with the LDA + U method by including a 
Hubbard- U term in the Hamiltonian for the Cd-4d and the 
Yb-4/ states. 
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